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1 Introduction 

The standard model neutrino interaction with other fermions is known to be extremely 
weak. Nevertheless, in some cases, even this weak interaction is crucial in the evolution of 
a neutrino system. One can recall the resonant amplification of neutrino flavor oscillations 
in background matter the Mikheyev-Smirnov-Wolfenstein (MSW) effect [1, 2] - which 
is the most plausible explanation of the solar neutrinos deficit [3]. Typically background 
matter, which a neutrino interacts with, is considered to either be at rest or move with 
a constant velocity. In both cases matter effects can influence the neutrino oscillations 
picture. 

Recently neutrino flavor oscillations in matter moving with an acceleration was dis¬ 
cussed in refs. [4-6]. The modification of the neutrino refraction index owing to a nonzero 
matter acceleration was accounted for in refs. [4, 5] without consideration of possible non- 
inertial effects. The neutrino motion in an accelerated frame accounting for noninertial 
effects was described in ref. [6]. The influence of the Unruh radiation on neutrino flavor 
oscillations was discussed in ref. [7]. 

Linearly accelerated particles reveal another interesting effect consisting in the emis¬ 
sion of a thermal radiation, which was first predicted by Unruh [8]. The temperature of this 
radiation appears to be proportional to the particle acceleration. This effect is still widely 
discussed in the literature (for a recent review see, e.g., ref. [9]). In ref. [10], it was shown 
that the Unruh effect strongly depends on the way to observe a thermal radiation. Never¬ 
theless there are some suggestions how to detect the Unruh radiation experimentally [11]. 
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In the present work we continue to study the influence of noninertial effects on elec- 
troweakly interacting particles. In refs. [ 6 , 12] we considered the electroweak interaction of 
neutrinos, electrons, and quarks with a rotating background matter in a corotating frame. 
In particular, the new galvano-rotational effect was predicted in ref. [12]. In this work 
we study the neutrino interaction with a linearly accelerated matter. Note that this kind 
of the matter acceleration can be implemented in various astrophysical media such as a 
supernova (SN) at the bounce stage [13] and jets from active galactic nuclei [14]. One can 
also expect intense neutrino beams in these environments. 

This work is organized as follows. In section 2 we recall the standard model neutrino in¬ 
teraction with background fermions which can move with a constant velocity. In section 3 
we derive and solve the Dirac equation, accounting for noninertial effects, for ultrarela- 
tivistic neutrinos interacting with a linearly accelerated matter. In section 4 we study the 
neutrino quantum states in a linearly accelerated matter. The neutrino-antineutrino 
pairs creation is described in section 5. In paticular, we reproduce the Unruh effect for 
neutrinos in electroweakly interacting matter moving with a linear acceleration and con¬ 
sider the correction to this effect owing to the specific neutrino interaction. In section 6 
we propose the astrophysical application of our results, consisting in the vv pairs creation 
in a SN explosion at the bounce stage. Finally, in section 7 we discuss our results. Some 
useful mathematical formulas are provided in appendix A. 

2 Neutrino interaction with matter in an inertial frame 

In this section we shall briefly remind how active neutrinos interact with background matter 
in the flat space-time. We describe the neutrino electroweak interaction in the Fermi 
approximation. 

The electroweak interaction of the flavor neutrino eigenstates a = e,/r, r, with a 
background matter consisting of electrons e, protons p, and neutrons n, can be described 
in the mean field approximation by the Langrangian 

-^int = ~ (2-1) 

a=e,/i,T 

where 7 ^ = 7 ^ (l — 7 ^) /2, 7 ^ = ( 7 *^, 7 ) are the Dirac matrices, and 7 ^ = i 7 '^ 7 ^ 7 ^ 7 ^. The 
effective current in eq. (2.1) has the form [15], 

J/t = V2GF + (2.2) 

f=e,p,n 

where Gp = 1.17 x 10“^ GeV~^ is the Fermi constant, is the hydrodynamic current of 
background fermions, Xj is the four vector of the mean polarization. The hydrodynamic 
current jj = UfU^ depends on the invariant number density Uf, i.e. the density in the 
rest frame of background fermions, and the macroscopic mean four velocity Uj of these 
fermions. The explicit form of in terms of the invariant density, four velocity, and the 

^ (1 2 ) 
invariant polarization (^f can be found, e.g., in ref. [16]. The coefficients in eq. (2.2) 
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(2.3) 


for Vf. have the form [15] 

„(/) _ r{/) _ C _1_ A _ A 

9l — -'lS ^y/4+0e/) 92 “ -'lS '^e/j 


( f) 

where I^s'^ is the third component of the weak isospin of type / fermions, Q/ is the value of 
their electric charge, ^ = sin^ ~ 0-23 is the Weinberg parameter, (5e/ = 1 for electrons 
and vanishes for protons and neutrons. To get </)• ’ for interactions with the same 
background fermions, we should set 6ef = 0 in eq. (2.3). 

The Lagrangian in eq. (2.1) can be obtained by the averaging of the fermionic vector 
and axial-vector currents in the rest frame of background fermions. Therefore, fermions, 
contributing to eq. ( 2 . 2 ), are supposed to have constant and homogeneous mean velocities. 
Otherwise (e.g, for matter moving with an acceleration; cf. section 3 below), eq. ( 2 . 2 ) is 
not valid. Indeed to obtain eq. (2.2) one should make a Lorentz boost to the rest frame of 
background fermions to make the averaging and then make a boost back to the laboratory 
frame. Such coordinate transformations are undefined if the Lorentz symmetry is broken. 

In our work we shall be mainly interested in the neutrino evolution in the dense electri¬ 
cally neutral nuclear matter of a protoneutron star (PNS). Let us assume that this matter 
is at rest and unpolarized. In this case, on the basis of eqs. (2.2) and (2.3) we get the 
effective potential as 




Gf 


(2.4) 


where Un is the neutron density. Note that the effective potentials for different neutrino 
types are approximately equal since the densities of charged particles inside PNS are much 
lower than that of neutrons. 

In the wake of the results of the recent experiments (see, e.g., ref. [17]), it is commonly 
believed that massive neutrino eigenstates, i.e. those which have definite masses, are the 
mixture of flavor neutrinos t'o defined above. It should be noted that the matter interaction 
of mass eigenstates is nondiagonal in the particle species. In the present work we shall not 
distinguish between mass and flavor eigenstates. For instance, this approximation is valid 
if we study the evolution of the t'e ~ W system, where the mixing angle @13 is known to be 
the smallest. Alternatively we may consider the neutrino system with the parameters far 
away from any MSW resonance. 

According to many theoretical models for the neutrino mass generation [18], the neu¬ 
trino mass eigenstates are likely to be Majorana particles. Nevertheless, in the present 
work, we shall assume that the neutrino mass eigenstates ipi, i = 1,2,..., are the Dirac 
particles. Note that, despite the great experimental efforts [19], the issue whether neutrinos 
are Dirac or Majorana particles still remains open. 


3 Dirac equation for a neutrino interacting with linearly accelerated mat¬ 
ter 

In this section we shall describe the motion of a neutrino in a linearly accelerated matter. 
Our analysis will be based on the exact solution of the Dirac equation accounting for 
noninertial effects. 
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The neutrino interaction with background matter described in section 2 implied the 
uniform matter motion with a constant velocity, i.e. in eq. ( 2 . 2 ) depends on neither 
coordinates nor time. The opposite situation would mean an accelerated matter. As was 
mentioned in section 2 , for accelerated matter we cannot assume a nonzero in 

eq. (2.1) and in the corresponding Dirac equation, written in the flat space-time, since the 
Lorentz symmetry is broken. 

Nevertheless, even for accelerated matter, we can always find a reference frame where 
matter is at rest. Assuming that matter is unpolarized, we get that only ^ 0 exists 
in this reference frame. Thus, in this frame, we can still use the effective potential defined 
in eq. (2.4) since it depends on the invariant density of background fermions. However, 
in this case the Dirac equation should be modihed to account for noninertial effects. This 
approach was developed in refs. [ 6 , 12 ], where we studied the propagation of electroweakly 
interacting particles (including neutrinos) in a rotating matter. 

It is known that the description of the particle motion by an accelerated observer is 
equivalent to the interaction of this particle with an effective gravitational field. Hence, in 
our treatment of the neutrino evolution, we should rewrite the Dirac equation in the curved 
space-time corresponding to a linearly accelerated frame, i.e in the Rindler space-time. If 
we suppose that the matter is accelerated along the z-axis, the interval in the Rindler 
space-time has the form [ 20 ], 

— dx^ — dy^ — dz^, (3.1) 

where is the metric tensor of the effective gravitational field and a is the proper accel¬ 
eration of matter. Here we use the Cartesian coordinates = {t,x,y,z). The interval in 
eq. (3.1) can be transformed to the Minkowskian one ds^ = dT^ — dA^ — dT^ — dZ‘^ by 
changing the coordinates 


T =zsinh at, 

1 , T 

t =—arctanh—, 
o Z 


Z =z coshat, 

>< 

II 

II 


Z =y/Z^-T^, 

II 

II 

(3.2) 


It should be noted that the metric in eq. (3.1) does not span over the all space-time. It 
covers only the sector Z > \T\ in the {Z,T) hyperplane. 

The Dirac equation for a neutrino interacting with background matter in a curved 
space-time can be obtained, using the results of refs. [ 6 , 12 , 21 ] and accounting for eq. ( 2 . 1 ), 
in the form, 

- m] = ^J^'jf.ix) [l - 7 ^(x)] (3.3) 

where is the neutrino bispinor, m is the neutrino mass, 7 ^(x) are the coordinate depen¬ 
dent Dirac matrices, = 9/^ + T^ is the covariant derivative, T^ is the spin connection, 
7 ®(x) = —(i/ 4 !)£;^'^“^ 7 ^(x) 7 jy(x) 7 Q(x) 7 ^(x), is the covariant anti¬ 

symmetric tensor in curved space-time, g = det{g^u), and is the effective external 
current of background fermions; cf. eq. (2.2). Note that, since we choose a comoving 
frame, only the zeroth component of is nonvanishing: = R 7 ^ 0, where V is given in 

eq. (2.4). In eq. (3.3) we omit the index Va for brevity: V = m = rriy^ etc. 
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One can check that the metric tensor in eq. (3.1) can be diagonalized using the following 
vierbein vectors: 


p M - 
Cn — 


(^^, 0 , 0 , 0 ^, = ( 0 , 1 , 0 , 0 ), e/ = ( 0 , 0 ,l, 0 ), = ( 0 , 0 , 0 , 1 ). 


(3.4) 


Indeed, the direct calculation using eq.(3.4) shows that rjab = where rjab = 

diag(l, —1, —1, —1) is the metric in a locally Minkowskian frame. 

To obtain the spin connection one should introduce the constant Dirac matrices in a 
locally Minkowskian frame as 7 “ = e“^ 7 ^(x). As shown in ref. [ 6 ], 7 ^(x) = i 7 *^ 7 ^ 727 ^ = 
7 ® does not depend on coordinates. The spin connection in the Dirac eq. (3.3) has the 
form [ 21 ], 


_ " ab _ i' 


(3.5) 


where Cab = (i/2) [ 7 a; 7 ?;]- are the generators of the Lorentz transformations in a locally 
Minkowskian frame and the semicolon stays for the covariant derivative. The explicit 
calculation on the basis of eq. (3.5) shows that the nonzero components of the connection 
one-form ujab = a^afe^dx^ are 


woiAt = -wiom = 0)- (3-6) 

Using eqs. (3.5) and (3.6) we get that i 7 ^(x)r^ = i 7 ^/ 2 z. 

Using the dehnition of 7 “, the Dirac eq. (3.3) takes the form, 

17 °^ + + i7^ ^02 - m 

Analogous Dirac equation in the Rindler wedge with U = 0 was studied in ref. [22]. Since 
eq. (3.7) does not explicitly depend on t, y, and z, we shall look for its solution in the form, 


-0 = -azj^V ( 1 - 7 '' 


'll). 


(3.7) 


if) = exp {-lEt + \pxX + ipyv) 'ipz, 


(3.8) 


where ipz = ipziz) is the wave function depending on z. 

It is convenient to rewrite eq. (3.7) as 

[7“Qa - m-h [/] ■i/’z = 0, (3.9) 

where = q°‘ — ^effAgg, ^eff is the effective electric charge, = {0,px,Py, —ifl^), 

(3.10) 


/ta _ 
^eff — 


1 / azV E i 


QeS 


is the potential of the effective electromagnetic field, and U = azV^^^^I2. 

Let us look for the solution of eq. (3.9) in the form, = [q^Qa + m — U]^, where <I> 
is the new spinor. The equation for 4> reads 


1 


E azV\‘ 


az 


■) 


-P± + 


a?‘z^V'^ 


— m 


/ E 

aV\ 

aV 


f E 

azV\ . 1 

+ ( — 2 

H - TT 

i«3- 

2z 

- - 



2 y 

2 


yaz 

2 J \ 


. 0„,5 


7 ° -|- mazV'j 7 ' 


$ = 0, (3.11) 
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where p]_= pi+ p^, as = and Sa = 7 ° 7 ^ 7 ®. 

The solution of eq. (3.11) can be found for ultrarelativistic particles. In the limit 
m —)• 0, we can represent <I> = in eq. (3.11), where (p = (p{z) is a scalar function and v 
is a constant spinor satisfying S 3 U = av and 7 ®u = xu, with a = ±1 and x = if) since 
both S 3 and 7 ® now commute with the operator of eq. (3.11). 

We shall be interested in the description of left active neutrinos satisfying (1 + 7^)^=0 
and having x = + 1 . One can show that, for right sterile neutrinos with x = “Ij the 
interaction with matter is washed out. Using the new variable p = we can write the 

equation for as 


P 


4 


where 


as p^\ + m? 

K = Ko + SKl-\—, Ko = ———, 

4 A\V\a 


K Pa = 0, 

(3.12) 

E 1 

(3.13) 



and s = sgn(U) = —1 for neutrinos in a neutron matter; cf. eq. (2.4). 


4 Classification of the neutrino quantum states 


In this section we shall find the solution of the wave equation for neutrinos interacting with 
accelerated matter in the explicit form and classify the neutrino quantum states. 

The solutions of eq. (3.12) have the form. 


CPo- 


Vp 


Vp 


(4.1) 


where and W\^^{z) are the Whittaker functions. Note that solutions of the Dirac 

equation, expressed via Whittaker functions, in the 1 + 1 de Sitter space-time were also 
obtained in ref. [23]. 

We shall suppose that the matter acceleration takes place in the spatial interval z\n < 
z < Zout, where Zin <C Zout- It corresponds to the change of the dimensionless variable 
Pin < P < Pont, where pin <C 1 and pout 1- Using eq. (A.l), we get that can 
be identified with the quantum state corresponding to a definite sign of the momentum 
projection along the z-axis at p ~ p-m, i.e. an in-state, whereas is the wave function of 
this quantum state at p ~ Pout- We shall call it an out-state. 

Therefore we can define two sets of the orthogonal neutrino wave functions: ±'il>n and 
(in and out solutions). These wave functions can be obtained from ^p and ±p using 
eq. (3.8) and the relation between and $ in section 3. We shall denote all the quantum 
numbers, like Px,y, E, and a, as n for brevity. 

The secondly quantized neutrino field 'll) can be decomposed using ±ipn and ^ipn in the 
form, 

^ [an(in)+V’n + &]^(in)_'0„ = ^ |^an(out)+V’n + blioutytpn , (4.2) 

n n 


where and bn are the annihilation operators of the quantum state corresponding to a 
definite sign of the momentum projection along the z-axis. These operators act on quantum 
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states in “in” and “out” Fock spaces. The vacuum in these Fock spaces is defined as 


dyi, (in, out) | 0 )in,out 


t)^d)|0)in,out — 0* 


(4.3) 


The creation and annihilation operators obey the usual anticommutation relations, 


On (in, out), a|j,(in, out) 


6 „(in,out), 6 )^, (in, out) 


= Sn 


(4.4) 


with all the other anticommutators being equal to zero. 

Now, let us find in tbe explicit form. For this purpose we represent the 

bispinor ■0^ of active ultrarelativistic neutrinos as ijjJ = (0,r/). Here we assume the Dirac 
matrices in the chiral representation [24] 


7 


0 



/ 0 cjfc \ 

( 0 ) ’ 



(4.5) 


where ak = (ci, < 72 , < 73 ) are the Pauli matrices. The two component spinors r^o- correspond¬ 
ing to the opposite helicity states with a = ± 1 , can be obtained if we choose the following 
constant bispinors Va'- 

uT = (l,0,0,0), uT = (0,l,0,0). (4.6) 

Finally, using eqs. (4.5) and (4.6), we get 773 - in the form. 


where 


v+ =n 




n =2y^^V\ 



E 

2a^ 


+ 1(7 


Vpdp + 



(cr^px) , 


(4.7) 


(4.8) 


is the projection operator, cr_i_ = (< 71 ,( 72 ), and p_L = {px,Py)- 

Let us first obtain out-states in the explicit form. Using eqs. (4.1), (4.7), and (4.8), as 
well as eq. (A.3), we can obtain the wave functions corresponding to the negative momen¬ 
tum along the z-axis. 


{ip) + /3+VF_3/4+iKo_iKi,-i/4-iKi(ip) 

\ -{Px+ ipy) lU-l/4+i<co-i«;i,l/4-iKi (ip) 

J_ ( -{Px- iPy) f^l/4+iK0-iKl,l/4-Hi<ci (ip) ) 

\/p y —iKi,l/4+iKi (ip) + /3-f^3/4-|-iKo-iKi -l/4+iKi(ip) J 


(4.9) 


Analogously one can derive the wave functions for the positive projection of the momentum 
as 


i^ I ®-|-iU./4—iKo+iKlil/4—iKi ( ip) A /^-|-i^3/4—iKo+i«;i,—1/4—iKi ( ip) | 

y/P\ -{Px + iPy) f^l/4-i<co-Hi<ci,l/4-iKi (-ip) J 

^ f {Px iPy) 1/4—iKo+iKi,l/4+iKi ( ip) 

y/P y 1/4—iKo+i«;i,l/4-|-iKi ( ip) T /3—3/4—iKo+irei,-l/4+iKi ( ip) 


(4.10) 
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The coefficients in eqs. (4.15) and (4.10) have the form, 


^0(^1 = — 2Ea 


ap 


'-/3± = 


=2v^e±3-/4 ^ i j , (4.11) 

where Q and C,' take the values ± independently. 

We define the norm of the total wave function as 

{'•PW) = j ^^x^/^PpE,p^,pyl^{x)1pE',p',,p'y = llV'f - E')5'^{-p± - p'±)> (4-12) 

where = ^^ 7 ^. The norm of a two component spinor 7 can be defined analogously. Using 
eqs. (4.9) and (4.10), we get that 


cji ( p-Y = + 7 +, 0-1 Yp-Y = p+. 

Therefore, one can establish the relation between the norms of the spinors, 


l|2 11+ ||2 

'n-W = w+ 


l|2 II- ||2 

= 'n+W ■ 


(4.13) 


(4.14) 


Note that the analog of eqs. (4.13) and (4.14) was found in ref. [25]. 

The in-states can be obtained in the same manner as eqs. (4.9) and (4.10). We just 
give the final result, 

_ ^ ( —®+-^— 1 / 4 +iKo—iKi,l/4—iKi (ip) T— 3/4+iKo~i^li~l/4—ifci (ip) j 

~ - (Px + iPy)M_i/4+i^Q_i^^^l/4_i«,^(i/9) )’ 


Vp\ 

1 f 


-V- = 


{Px iPy)-^l/4+i«;o~i^lil/4+i«:i (ip) 


\/P y —-^l/4+i«;o~i^ld/4+iKi (ip) T —■^3/4+i«;o~i^lj~l/4+iKi (ip) j 


and 


(4.15) 


_ ( +^+^1/A—iK,Q+iKi,l/A— ikY ip) T A—ikq+iki,—1/A—ikY ip) 

- (Px+ ipj/)Mi/4_i^Q+i^^^i/4_i^^(-ip) 


+P- = 


Yp 

u 


{Px iPy) 1/4—iKo+i^i;l/4+i/^i ( ^P) 


y/P Y+^—1/4—iKo+i«^i,l/4+i«:i ( ^p'} +/5—3/4—i«:o+i^i!~l/4+iKi ( ^p) 

The coefficients in eqs. (4.15) and (4.16) are 


(4.16) 


^a^,=-2EJ^P, _/3± = (4«^i ± i) , 

\ ap 


-I3± =e--/^v^(4^i±i), 


(4.17) 


where C and (' again take the values ± independently. Analogously to the out-states, we 
can establish the following relations between the norms of the wave functions: ||_ 7 _|p = 








Il+r^+ll^ and ||+77_||^ = ||_7/+|p. These expressions are valid since it follows from eqs. (4.15) 
and (4.16) that ai (-?]-)* = +t/+ and ai (+t/_)* = _r/+. 

To complete the determination of the states we need to establish the relation between 
the norms of the wave functions for the same helicity and different signs of the momentum 
projection. This relation can be found in the case when ki 3> kq- Using eqs. (A.4) 
and (A.6), we get the ratio of the norms of and as 


7?| 


vre 


—27r/^i 


-77_| 


cosh 27rKi 


(4.18) 


Note that one can derive the analogous relation between, e.g., the wave functions +? 7 _ and 
_r/_. However, we shall not present it here. 


5 Particle creation in an accelerated matter 

In this section we describe the creation of pairs in a background matter moving with a 
linear acceleration. We also reproduce the Unruh effect. 

Both out- and in-states are the complete sets of orthogonal wave functions. Thus we 
can expand any in-state wave function over the system of out-states, 

= '^ G C = ±, (5.1) 

C'=± 

where the Bogoliubov coefficients G satisfy the unitarity conditions, 

E o (''!<) G (cl^") = E G (<'t) G («k") = in- (5.2) 

C=± C=± 

Here g(cK') = are the coefficients of the inverse transformation. 

The relation between in and out wave finctions in eq. (5.1) is equivalent to the Bogoli¬ 
ubov transformation of the creation and annihilation operators defined in eqs. (4.2)-(4.4), 

an(out) =G (+1+) an(in) G (+|_) 6j,(in), 

^n(out) =G ("1+) a„(in) G ("I-) 6]j(in). (5.3) 

Therefore, the nonzero mean value of the number density operator of out-states over the 
“in” vacuum, in(0|an(out)a„(out)|0)in = |G('''|_)|^ / 0, can be interpreted as the creation 
of a lyp pair. 

From the formal point of view the vacuum instability leading to the pairs creation, 
described in the present work, is analogous to the pairs creation by an external field de¬ 
pending on spatial coordinates. It should be noted that a more detailed analysis of this 
problem was recently presented in refs. [26, 27]. 

Setting C = — in eq. (5.1) we get that 

=G'(+|_)+r/^-hG("|_)"r/^. (5.4) 
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The Bogoliubov coefficients in eq. (5.4) are related to the scalar products, defined in 
eq. (4.12), of the corresponding in- and out-states as 

, G{~\-) = . (5.5) 


Since we study ultrarelativistic particles, there is a strong correlation between their mo¬ 
mentum and spin: neutrinos are left-polarized, whereas antineutrinos are right-polarized. 
If cj = — 1, the nonzero coefficient G ("’"I-) corresponds to the creation of the pair of v with 
{pz > 0, a = —1) and P with {pz < 0,a = —1). In case a = +1, G('''|_) describes the 
creation of v with {pz <0,a = -|-1) and u with {pz > 0,a = -|-1). However, these two sets 
of states are identical. Thus, if we define the corresponding probabilities 

P. = |G(+|_)|^ (5.6) 


we should get that P- + = 1. In fact, this relation between the probabilities is the 

consequence of eq. (4.13). 

To derive in the explicit form, we should use eq. (A.5). For cr = —1, on the basis 
of eqs. (4.9), (4.10), and (4.15), we have 


_r/_ =r ( - -h 2iKi ) exp 


ITT 


7r(Ko - «^i) - ^ 


exp (—TTKi -|- 37ri/4) _ 

T (1 iKo) 

Substituting eq. (5.7) to eq. (5.5) we obtain 


P- + 


1 


T ( 1/2 -|- 2 i«;i — iKo) 


p- 


G( |_) / 37ri\ II u_|| T {1/2-\-2\ki — ikq) 

- = -P + l-j ||.,_||^ r(i + i.o) ■ 


G(+l-) 

Analogously, if a = -|-1, one gets that 


/3 \ 

m 

1 - - 2i/ti 1 exp 

7r{Ko - Ki) + — 


X 


exp (vTKi -|- 37ri/4) _ 

T ( 1/2 — 2 ifi:i -|- iKo) 


P+ + 


T (1 - i/to) 


P+ 


and 


G'( I-) f 37ri\ II p+\ 

= exp vr/ti -|-— 


r (1 - i/to) 


G(+|_) 4 ; ||+^^fr(l/2-2i/ti+i/to)' 

Finally, accounting for eqs. (5.2), (5.6), and (A.6), we get that 


P_ = 


P^ = 


1 + e 


—27rKi II I 


sinh vr/to 


|+r/_|p /to cosh7r(2/ti — /to) 


-1 


1 + e 


2nKi II P+W Kocosh7r(2/ti - /to) 


1 “1 


'p+\ 


sinh TT/to 


(5.7) 


(5.8) 


(5.9) 


(5.10) 


(5.11) 


Taking into account eq. (4.14), one can see that P- + = 1, as it should be. 
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Let us discuss the situation when ki S> 1 S> kq, i.e. the creation of vv pairs is mainly 
owing to the noninertial effects. We can define the number of created pairs per second in 
the energy interval dE as dA^ = P-dE/2 ti. On the basis of eqs. (4.18) and (5.11), we get 
the distribution of created particles in this limit as 


diV 

dE 


= exp 


E + 5E\ 

Teff ) ’ 


?eff 


a 


5E 


p\ + m? 

8 | 4 ^| 


(5.12) 


Here we reproduce the temperature of the Unruh radiation Teg [8]. The correction to the 
Unruh effect owing to the specific electroweak neutrino radiation is given by the quantity 
5E in eq. (5.12). Note that eq. (5.12) is obtained under the assumption that E S> dE. 


6 Neutrino pairs creation in SN 


In this section we shall discuss a possible application of the obtained results for the de¬ 
scription of the ui) pairs creation in a core collapsing SN. In particular, we shall obtain the 
astrophysical upper limit on the neutrino mass. 

The distribution of created uv pairs is given in eq. (5.12). If we discuss the case of 
the neutrinos propagation along the z-axis, the correction to the Unruh effect does not 
suppress the pair creation if 6E <C This condition can be rewritten as the constraint 
on the neutrino mass: m Trier, where 


rric 


= 2 


|U|a 


TT 


( 6 . 1 ) 


Note that the constraint in eq. (6.1) is theoretical. It means that, if vi? pairs, created 
because of the vacuum instability in an accelerated matter, are observed experimentally, 
then the neutrino mass should be less than nicr in eq. (6.1). 

Let us evaluate rricr for a core collapsing SN. When the plasma pressure cannot equili¬ 
brate the gravity of the upper layers of a protostar, the radius of a protostar core starts to 
decrease until the matter density reaches pn ~ 10^^ g • cm“^. Assuming that PNS mainly 
consists of neutrons, this value corresponds to the neutron density ~ 6 x 10^^ cm~^. At 
this point no further core compression happens and initially falling matter bounces from a 
dense core which then turns out to be PNS. According to the results of ref. [13], the matter 
velocity changes by lAuj ~ 5 x 10® cm • s~^ during At ~ 1ms after a bounce. It gives 
the nonzero matter acceleration a = lAuj/Af ~ 5 x 10^® cm • s“®. Using eq. (2.4) and the 
chosen values of all the parameters, we get that rricr ~ 7.2 x 10“^ eV. The obtained value 
of rricr is comparable with the (theoretical) upper bounds on the neutrino masses derived 
in ref. [28]. 

The total number of vu pairs emitted during At from the whole surface of PNS can 
be obtained by integrating eq. (5.12) over E and p_L. The final result reads 

" -i 


N = 8\V\T^sAtR^ exp - 


m 


S\V\TcsJ ’ 


( 6 . 2 ) 


where R is the PNS radius. The maximal matter acceleration at a bounce happens at 
R ~ 10^ km [13]. Therefore, assuming that m ~ rricr, on the basis of eq. (6.2) one gets that 
N ~ 10^^ neutrino pairs can be emitted by PNS. 
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This vv flux is much less than the flux of Vf. emitted in a neutrino burst, which can be 
estimated as ~ 10^^ [29] . However, the typical energy of v and i>, created owing to the 
vacuum instability in a dense nuclear matter in PNS, is in the eV range; cf. ref. [28], and 
lOMeV. Therefore, at the bounce stage, high energy i/g’s will be absorbed by the 
dense matter of PNS whereas low energy v and v will freely escape PNS. Indeed, using the 
neutrino scattering cross sections given in ref. [30], one finds that the mean free path of v 
and V, with energies in eV range, in background matter with the density pn ~ 10^^g-cm“^ 
is about 10^® km, which is far beyond R = 10^ km. Thus one can consider the emission of 
vv pairs created by the accelerated matter in PNS as a precursor of the neutrino burst. 

Nevertheless, it is unlikely that such v and P, created by the proposed mechanism, 
can be detected using modern experimental techniques. Despite the recent suggestions to 
observe nonrelativistic relic neutrinos with sub-eV energies [31], the flux of vv pairs emitted 
in a SN explosion is quite low to be detected in any existing neutrino telescope even if this 
explosion happens in our Galaxy. 

7 Conclusion 

In conclusion we mention that, in the present work, we have studied the evolution of 
neutrinos electroweakly interacting with a linearly accelerated background matter. We 
have briefly reminded the standard model neutrino interaction with background fermions 
in the flat space-time in section 2. Then, in section 3, we have derived the Dirac equation 
for a massive neutrino, electroweakly interacting with background fermions, in the Rindler 
space-time corresponding to a linearly accelerated matter. The solution of this Dirac 
equation has been found for ultrarelativistic particles. The quantum states of neutrinos 
have been analyzed in section 4. In section 5, we have applied our results for the description 
of the vv pairs creation in an accelerated matter. Finally, in section 6, we have considered 
the vv pairs creation in a collapsing SN at the bounce stage. 

As a main tool in our study we have used the method of exact solutions of the Dirac 
equation for a massive neutrino rewritten in a comoving frame, where matter is at rest. 
Firstly, this method enables one to unambiguously determine the neutrino interaction 
potentials with background fermions which move with an acceleration [6, 12]. Secondly, 
unlike refs. [4, 5], where the matter acceleration was accounted for only in the neutrino 
refraction index, our approach allowed us to examine noninertial effects. 

We have demonstrated that the neutrino interaction with a linearly accelerated matter 
can induce the vacuum instability which leads to the creation of vv pairs. We have discussed 
the situation when an observer was in the accelerated frame since we used the coordinates 
in the Rindler wedge. Therefore, the main contribution to the vv pairs creation is owing to 
the Unruh effect [8]. We have rederived the temperature of the Unruh thermal radiation 
and studied the correction to the Unruh effect due to the specific electroweak neutrino 
interaction. 

The electroweak decays of accelerated particles were studied refs. [32, 33]. The prob¬ 
ability of the decay of an accelerated proton p e^nv^. was shown in refs. [32, 33] to 
be nonzero. In those works the hadronic current was described quasiclassically in the flat 
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space-time. Moreover it was demonstrated that the decay rate is equal to the sum of the de¬ 
cay rates of cross symmetric processes. These cross symmetric processes were studied using 
the exact solutions of the Dirac equation in the Rindler wedge. The results of refs. [32, 33] 
may be interpreted as the indirect proof of the existence of the Unruh effect. 

Unlike refs. [32, 33], in the present work we consider the elastic forward neutrino 
scattering off accelerated background fermions / with a possible creation of vv pairs; vf A 
We treat background fermions as a classical external current. Therefore we confirm 
the result of refs. [32, 33] that the Unruh effect for neutrinos takes place if background 
fermions are treated quasiclassically. 

In our analysis we use the exact solution of the Dirac equation which accounts for both 
the interaction with matter and its acceleration, i.e. our results exactly take into account 
all terms in the perturbative expansion over Gp. It is the main difference from refs. [32, 33], 
where the electroweak decays were studied on the tree level linear in Gp. The chosen way 
to describe neutrinos interacting with an accelerated matter allowed us not only establish 
the existence of the Unruh effect for these particles but also to obtain the correction to 
this effect in eq. (5.12). 

The obtained electroweak correction to the Unruh effect was shown to contribute to the 
uv pairs creation in a core collapsing SN. Assuming that the process of the pairs creation 
is not suppressed, we have obtained the upper (theoretical) limit to the neutrino mass; 
cf. eq. (6.1). This constraint is in agreement with the result of ref. [28]. Unfortunately, 
nowadays it is almost impossible to experimentally detect u and v created in the proposed 
mechanism. 

It should be noted that previously the creation of lyv pairs in dense matter of a neutron 
star (NS) was considered in refs. [34-36] using the Schwinger mechanism. For the first time 
the study of this problem based of the quantum field theory was undertaken in ref. [28]. 
However, in ref. [28] the case of the time dependent effective potential of the neutrino 
interaction with the background matter was discussed. Thus the model considered in 
ref. [28] corresponds to the neutronization stage of PNS. 

The process of the uv pairs creation described in the present work can be attributed 
to the dependence of the external field (background matter) on spatial coordinates. The 
general description of the vacuum instability in case of coordinate dependent external 
fields was recently discussed in refs. [26, 27]. In our work for the first time we consider 
the coordinate dependent effective potential of the neutrino interaction with background 
matter as a source of the neutrino vacuum instability. 

Finally we mention that the creation of vv pairs considered in the present paper is 
different from the bremsstrahlung emission of vi> pairs in the NS matter; cf. refs. [29, 37]. 
Neutrino pairs emitted in the latter effect have energies proportional to the temperature 
of the nuclear matter in NS and thus correspond to high energies. On the contrary, we 
predict the coherent emission of low energy particles. 
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A Properties wi v* 
In 


this appendix we list some of the useful properties of Whittaker and Euler gamma 
Lctions. These expressions are taken from ref. [38]. 


functions. These expressions are taken from ref. 

The asymptotic expansions of Whittaker functions read 

-h ..., at l^l -)• 0 , 
Wx,p.{z) ~ -k ..., at |z| —)• 00 . 

The derivatives of the Whittaker functions have the form, 


dW±i^^^{±ip) e^3iW4 


dp 

Vp L' 


/2p- 


V 2p 

dM±iK,^(±ip) 


dp 

Vp 

dIU±iK,^(±ip) 

g±i7r/4 

dp 

Vp 

dM±iK,^(±ip) 



«^=Fi X 


^±iK-l/2,At-l/2(=tip) 


K-l/2,M-l/2(±ip) - ^ ± 2 ) 

■^±i«;+l/2,M-l/2(±ip) - ^ ^2p ^ 2 ) 


dp ^/p 

At p = A ± 1 / 2 , Wx^^{z) takes the form, 

W"a,a-i/ 2 (^) = Wx,x+i/ 2 {z) = z-^e^/^T{l + 2 A, z), 

where T{a,z) is the incomplete gamma function. 

The function Mx,^{z) is related to Wx,^{z) by the following expression 

W x,^.{e-^ 


Mx,im{z) = r (2p -k 1) e 


—IttA 


r (1/2 + // — A) 


I i7r(At+l/2) jz) 

r(i /2 + p + A) 


(A.l) 

(A.2) 


(A.3) 


(A.4) 


(A.5) 


which is valid if — 7 r /2 < argz < 3 ' 7 r /2 and 2p 7 ^ —1, —2,.... 

The absolute value of the Euler gamma function of the complex argument in some 
cases can be evaluated as 


r 



2 


TT 

cosh TTX ’ 


r(i + ix)|2 


TTX 

sinh TTX ’ 


(A.6) 


where x is a real parameter. 
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